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ABSTRACT 

We explore the possibility of resolving an image of a damped Lya (DLA) system in absorption against an 
extended, diffuse background X-ray source. Typical columns of neutral hydrogen in DLA systems are high 
enough to block out up to ^ 30% of the soft X-ray flux at an observed photon energy of 0.5 keV, and we find 
that ~ 1% of the area of extended X-ray sources at z > 1 have their 0.5 keV flux reduced by at least 20% because 
of intervening DLA systems. We discuss the observability of such absorption and find that > 300 photons per 
angular resolution element are required in the 0.3-8 keV band for its detection, and in order to distinguish it from 
intrinsic surface brightness fluctuations. For the surface brightness of the currently known high-redshift extended 
X-ray sources, this requires an integration time of a few Msec on Chandra if the maps are smoothed spatially to 
w 2" resolution. The exact required integration time depends on the DLA system's column density, metallicity and 
most strongly, its redshift. Current X-ray telescopes are likely to detect DLA systems with A^hi < 10^^ cm"^ only 
out to z w 2.3. The availability of DLA systems with a suitably high column-density for a silhouette detection is 
currently poorly known. We suggest that at low redshifts, archival data of bright X-ray point sources may be useful 
in constraining the high-A^ni end of the column density distribution. We briefly discuss an alternative strategy of 
searching for extended X-ray sources behind known DLA systems. Although with current X-ray telescopes the 
detections are challenging, they will be within the reach of a routine observation with a next generation X-ray 
telescope such as the X-Ray Evolving Universe Spectrometer (XEUS) or Generation-X, and will deliver novel 
constraints on the nature of proto-galaxies. 

Subject headings: cosmology: theory - galaxies: formation - galaxies: high-redshift - quasars: absorption lines 
- X-rays: galaxies 



L INTRODUCTION 

Damped Lya (DLA) absorbers owe their name to the pres- 
ence of broad damping wings in their absorption profiles caused 
by column densities of neutral hydrogen in the range A^ei ~ 10^" 
to 5 X 10^'cm"^. These high column densities are only found in 
the cold disks of late-type galaxies in the local universe, but to 
identify the galaxies hosting this gas at cosmological redshifts 
has proven difficult, because of the overwhelming brightness of 
the background quasar in the optical. Only a dozen probable 
hosts for DLA systems at z < 1 have been identified and their 
morphologies span a wide range (Chen & Lanzetta 2003; Rao 
et al. 2003). 

So far several hundred DLA systems have been found in op- 
tical and ultraviolet quasar absorption spectra in the redshift 
range z ^ 0.1-4.6 (Curran et al. 2002). It has been shown that 
at z > 3 they are the main contributors to f^Hi (Rao & Turnshek 
2000; Storrie-Lombardi & Wolfe 2000; Prochaska & Herbert- 
Fort 2004). and make up the largest reservoir of cold gas avail- 
able for starformation. Indeed, ilni at z ~ 3 is comparable to 57* 
at z ~ 0, the cosmic mass density of stars in the local universe 
(Wolfe et al. 1995). 

It has therefore been argued that DLA systems are the pro- 
genitors of present-day galaxies. The exact nature of these pro- 
genitors however, is controversial: some have argued that DLA 
systems are large proto-galactic disks, based on (1) the asym- 
metric profiles of the low-ionization absorption lines, which are 
explained if the metals are present in a large, rapidly rotating 
disk (Prochaska & Wolfe 1997) and (2) the large (R > 15-30 
kpc;-e.g., Chen & Lanzetta 2003, Chen 2005) physical sepa- 
ration between the absorber and possible associated starlight 
in a few cases. Others have argued that these data are consis- 
tent with DLA systems being a collection of much smaller (< 5 
kpc) compact HI clouds falling onto a galaxy in their vicinity 



(Mailer et al. 2001;Haehnelt et al. 1998) and that these models 
provide a better fit to the kinematic data on the high ionization 
gas (Mailer et al. 2003). The collection of clumps in this sce- 
nario extends over an angular scale comparable to that of the 
single large disk. A distinguishing feature between the models 
is therefore the smoothness of the DLA system's gas. 

Determining the nature of the progenitors of our present-day 
galaxies is an important step in understanding galaxy forma- 
tion in general. Obtaining a spatially resolved image of the 
neutral gas in a DLA system will be an important step in set- 
tling this issue. Unfortunately, detecting the neutral hydrogen 
in DLA systems in emission at z > 0.5 using the 21cm line 
on an existing telescope such as the Giant Metre-wave Radio 
Telescope (GMRT) is not feasible (e.g. Bagla 1999). Among 
planned future radio telescopes, only the Square Kilometer Ar- 
ray (SKA) has the appropriate sensitivity and frequency range 
to make such observations. 

A major hindrance for constraining the properties of the DLA 
system's gas is the point-like character of the background quasars. 
Finding DLA systems in absorption against extended background 
sources would provide valuable additional information. In one 
example, Briggs et al. (1989) found that the dimensions of a 
disk-like absorber at z = 2 are > 1 1 kpc against an extended 
background radio source; ruling out compact gas configurations 
such as dwarf galaxies. In another case, Kanekar & Briggs 
(2003) have found HI 21cm absorption in a gravitational lens 
system at z = 0.76 toward a quasar, from which the column den- 
sity of hydrogen in the lens was estimated to be (2.6 ± 0.3) x 
10^' cm"^. The lens is therefore a DLA system, but because the 
largest separation between a pair of radio images was small, 
~ 0.7" , the gas configuration and extent could not be con- 
strained. Fynbo et al. (1999) detected extended Lya emission 
from a DLA system, that is illuminated by two nearby (Az ^ 
0.01) quasars, and found the emission to extend over a region of 
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3" X 6". A disadvantage of this technique is that it only works 
when a bright source of ionizing radiation is physically very 
close to the DLA system. 

In this paper, we examine the complementary possibility of 
detecting DLA systems against an extended background source 
emitting X-rays. The motivation for this approach is the dis- 
covery of the extended X-ray source 4C 41.17 by Scharf et al. 
(2003) at z = 3.8, with diffuse emission extending over ^ 100 
arcsec^. The X-ray spectrum is non-thermal and has a power- 
law index consistent with that of radio synchrotron emission. 
Scharf et al. (2003) conclude that the X-rays are most likely in- 
verse Compton scattered cosmic microwave background(CMB) 
or far-infrared background photons. A similar source, 3C 294, 
was found by Fabian et al. (2003) at z = 1.8 with diffuse emis- 
sion extending over ~ 200 arcsec^. 

A silhouette of a galaxy containing neutral hydrogen against 
a background X-ray cluster has akeady been found at low red- 
shift (z = 0.02, Clarke et al. 2004) and encourages us to explore 
diffuse X-ray sources as possible background screens against 
which to image DLA systems. Bechtold et al. (2001) found 
strong evidence for X-ray absorption in the direction of a z = 0.3 
DLA system in the spectrum of a z = 1.1 quasar, a bright X-ray 
point source with one faint jet (Siemiginowska et al. 2002). 

The outline of the rest of this paper is as follows. In §2, 
we introduce the X-ray optical depth through a DLA system 
and calculate the fraction of the sky covered by DLA systems 
as a function of their value of this optical depth. In §3, we 
describe how absorption by a DLA system can theoretically be 
distinguished from intrinsic surface brightness fluctuations and 
calculate the required integration times on Chandra to do this 
observationally. In §4, we discuss various contaminants that 
could mimick absorption by a DLA system; we also address 
the most important model uncertainties. The "inverse method" 
of searching for X-ray emission behind known DLA systems is 
discussed in §5. Finally, in §6, we present our conclusions and 
summarize the implications of this work. 

Throughout this paper, we focus on an observed photon en- 
ergy of £ = 0.5 keV. This choice reflects the fact that the absorp- 
tion signature is strongest in the range E ^ 0.3-0.6 keV (see 
§3), and that X-ray detectors are significantly less efficient at 
lower energies. We adopt the background cosmological param- 
eters n,„ = 0.27, Ha = 0.73, and h = 0.71 (Spergel et al. 2003). 

2. SILHOUETTES OF DLA SYSTEMS 

2.1. The X-Ray Optical Depth ro.5. 

Typical neutral hydrogen column densities in DLA systems 
are in the range A^hi ~ 10^** to 5 x 10^' cm"^ and are large 
enough to block a non-negligible fraction of X-rays emitted 
by sources that lie behind them. The optical depth through a 
column A^Hi of hydrogen atoms for a photon of energy E is 
t{E) ~ 0.43 (A^Hi/10^^)(£'/0.5 keV)-^^ ^, for a gas that consists 
solely of a cosmic mix of 76% neutral hydrogen and 24% he- 
lium by mass. Observations indicate that the mean metallic- 
ity of DLA systems is ^ 10% of the solar value with a mod- 
est redshift evolution (Prochaska et al. 2003). The presence of 
metals increases the optical depth. For example, the cross sec- 
tion per hydrogen atom in a gas with metallicity Z = O.IZq is 
~ 1.2(£'/0.5 keV)" -' times larger than in a metal free gas (e.g. 
Wilms et al. 2000). This shows that DLAs can show up in ab- 
sorption against bright background X-ray sources (as detected 
already for a point-source by Bechtold et al. 2001). 

We define the parameter T0.5 as the optical depth through 



the DLA system at the photon energy E = 0.5 keV in the ob- 
server's frame, ro.5 = A^HiO'tot(0.5(l +z) keV), where z and A^hi 
are the redshift and column density of neutral hydrogen atoms 
in the DLA system, respectively, and crtot(£) is the total absorp- 
tion cross-section per hydrogen atom. This cross-section in- 
cludes absorption by helium and heavier elements. To calculate 
crtot(E) of the DLA gas, we use the model described by Wilms 
et al. (2000). This model includes up-to-date photoionization 
cross sections and heavy elements in a solar abundance pattern, 
as well as H2 molecules. We note that X-ray absorption arises 
predominantly from inner shells ionizations of metals, and ion- 
ization corrections are only important at gas temperatures much 
higher than those expected to occur in DLA systems. 

2.2. The Influence of DLA Metallicity on atot and T0.5. 

We consider the specific class of DLA systems with T0.5 = 
0.1. Given crtot(E) and requiring A^HiO'tot(£' = 0.5(1 +z) keV) = 
T() 5, we obtain a column density of hydrogen Nmiz) for these 
DLA systems. In Figure 1, we show Nyaiz) as a function of 
redshift. 
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Fig. 1 . — Column density of hydrogen required to produce an absorption 
optical depth of ro 5 = 0.1 for photons with observed energy E = 0.5 keV, as 
a function of redshift for three metallicities. To illustrate the importance of 
metallicity in DLA systems, we show three cases: DLA systems with no metals 
(dashed-blue line), with Z = Zq (dotted-red line) and Z = 0. IZ0 (solid-black 
line). The last is our fiducial model. The highest observed column of neutral 
hydrogen among known DLA systems is log(A'Hi/cm"^)=21.7 (Curran et al. 
2002). 

To illustrate the importance of the assumed metallicity, we show 
three cases: Z = (dashed-blue line), Z = Zq (dotted-red line 
), and Z = O.IZ0 (solid-black line). The last is our fiducial 
model, based on the mean observed metallicity. Note that the 
column of hydrogen required to produce a certain value of tq g 
can vary by an order of magnitude depending on the metallic- 
ity of the gas. The reason A^hi drops slightly but suddenly at 
z ^ 0.4,0.7, 1.6,... is that E = 0.5(1 +z) keV at these redshifts 
equals the K-shell ionization energy of one of the heavy ele- 
ments and cTtot 'jumps' to a larger value. For any other value of 
ro.5, the corresponding column can be obtained by multiplying 
the value taken from Figure 1 by ro.5/0.1. We find that the fol- 
lowing formula provides an accurate fit to the relation between 
A^Hi, Z, Z and ro.5 for a DLA system: 
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To.5 = 0.014 



A(Z) ( 



1021 cm- 
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-3.2+B(Z) 
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\0.91n 



where A(Z) = 1 .0+5.39Z, and B(Z) = 0.64 arctan [(Z/0. 19) 

This formula is good to within 5% for Z < O.IZ0, and to 
within 10% for Z < Zq. The main reason for the deteriorating 
accuracy of the fitting formula at larger metallicities, is that the 
'jumps', in the absorption cross-section are more pronounced 
in these cases. It is exactly at these jumps that the fit is worse. 

2.3. The Covering Factor of DLAs. 

In §2.2 we calculated the column of hydrogen in a DLA sys- 
tem required to produce a given value of the X-ray optical depth 
To, 5 as a function of its redshift. In this section we calculate 
what fraction of the sky is covered by DLA systems with vari- 
ous minimum values of T0.5. 

The total number of DLA sytems along a single line of sight 
(LOS) with To 5 or larger between redshift and Zs, where Zs 
is the redshift of a hypothetical extended background source, is 
given by: 



Af= / dz 



dNuifiNHuz), 
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where Nhi(z) is the column of hydrogen corresponding to ro,5 
and /{Nhi , z)dNmdz is the number of DLA systems in the range 
A^Hi±'^A^Hi/2 and z±dz/2. We adopt /(Nhuz) from Peroux 
et al. (2003, hereafter P03), who find that the data can be best 
fitted with a F function used in studies of the galaxy luminos- 
ity function.' Note that JV can equivalently be interpreted as 
the two-dimensional covering factor of these DLA systems; 
/cov = A/". The exact form of f{NHi,z) is an important uncer- 
tainty in our model. The total sample of DLAs is still too small 
to constrain /(A^hi,z) well, especially at the low redshifts and 
high-column densities that are most important for our purposes. 
We discuss this uncertainty in more detail in §4.2 below. 

In Figure 2, we plot /cov as a function of ro.5, with the X- 
ray source located at = 3.8 (motivated by the source found by 
Scharf et al. 2003, 4C 41.17). The three thin curves represent 
the same three cases shown in Figure 1 in which the DLA gas 
has no metals (dashed-blue line), Z = O.IZq (solid-black line) 
and Z = Zq (dotted-red line). The covering factor increases 
with metallicity because the column of hydrogen required to 
produce a certain value of T0.5 decreases with increasing metal- 
licity. (The thick solid line addresses a scatter in metalhcities, 
as explained in §2.4 below.) 

Figure 2 shows for example, that for the fiducial model (Z = 
0. IZq) a fraction /cov ~ 0.03 and 0.01 of the sky is covered by 
DLA systems that absorb at least ^ 10% and ^ 20%, respec- 
tively, of the 0.5 keV photons from the background object. For 
different energies, te ^ T0.5 x [(Ttot(£')/o'tot(0.5 keV)], since the 
absorption cross-section is close to a pure power-law. 

The main conclusion from Figure 2 is that for a source sim- 
ilar to 4C 41.17, which has a size of ~ 100 arcsec^, observed 
at ^arcsecond resolution, one can hope to find several pixels 
covered by DLA systems with opacities up to tq s ^0.1 -0.2. 
Of course, for a source covering a larger solid angle and/or in 

'Note that we obtain f{Nm,z) by multiplying the functions quoted in P03 by 
dX{z)/dz- Their function describes the number of DLA systems in the range 
X{z)±dX{z}/2, where dX(z)/dz = {l+z)^iHo/H{z)) is the absorption distance, 
originally introduced by Bahcall & Peebles (1969), where Hq and H(z.) are the 
Hubble parameters at redshift and z, respectively. 



the case of an unusually metal-enriched DLA system, there is a 
significant probability of finding more opaque absorption. For 
example, in the case of an X-ray emitting galaxy cluster with ^ 
10 times the area of 4C 41 . 17 (see below), the figure shows that 
DLA systems with ro s > 0.3 should typically obscure several 
pixels. 
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Fig . 2. — Covering factor /cov on the sky of DLA systems between z = and 
Zi = 3.8 that have an X-ray optical depth tq 5 or larger, plotted as a function of 
TO 5 . The three curves represent the same three metallicities as in Figure 1 . The 
thick, solid line represents an alternative model, discussed in §2.4, in which the 
DLA systems have a log-normal distribution with mean log(Z/ZQ ) = — 1 .0 and 
o-iogZ/ZQ = 0.5. 



2.4. Uncertainties in /cov 

The calculation of /cov has two main uncertainties. They are 
discussed here. 

The larger of the two uncertainties is associated with the 
function f{Nui,z)- The redshift of the source, Zs, has only a 
small effect on the results for To 5 > 0.1. The reason is that 
with the f(Nni,z) adopted from P03, the absorption is domi- 
nated by low redshift systems; e.g. 75% of the contribution 
to /cov comes from z < 0.61 and the contribution from sources 
with z > 1.4 is practically negligible. This inference, however, 
depends sensitively on the poorly known redshift-evolution of 
f(Nm,z), and the contribution to /cov from high-z systems may 
be much larger (discussed further in §4.2). 

The second uncertainty involves the metallicity of the DLA 
gas. Figure 2 shows that the covering factor depends quite 
strongly on the assumed metallicity. Our fiducial model as- 
sumes Z = O.IZq, based on the mean observed value. There is 
however, considerable scatter in Z around this mean; the metal- 
licity can easily vary by an order of magnitude between individ- 
ual objects (Prochaska et al. 2003). Since enhanced metallicity 
increases /cov more than reduced metallicity decreases it, one 
may expect that incorporating the scatter enhances /cov. 

We quantified in an alternative model the effect of this ob- 
served scatter in metallicities on /ov A log-normal distribu- 
tion of metallicities around the mean log(ZDLA/ZQ) = -1 .0 with 
a standard deviation of criog(ZDLA/ZQ) = 0.5 was assumed (we im- 
posed a maximum value of Z = Zq since larger metallicities 
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have never been observed Prochaska et al. 2003). In this model 
< 10% of the DLA systems have metallicities log(ZDLA/Z0) < 
-1.6 and > -0.4. Our results for this alternative model are 
shown by the thick solid line in Figure 2. We find that the scat- 
ter causes a net small increase in /cov, as expected, but that in 
the entire range in ro.s plotted, the increase is < 25%, which is 
well within the other model uncertainty mentioned above. 

3. OBSERVABILITY 

The ideal background source against which to detect a DLA 
system is a bright and perfectly smooth extended X-ray source. 
In this idealized case a foreground DLA system (or that part of 
the DLA that intersects the extended X-ray source) shows up in 
absorption against the background source in the very soft X-ray 
band (0.3-0.6 keV). In practice, extended X-ray sources are not 
smooth, but have surface brightness fluctuations that are much 
larger than the ^ 10% level introduced by absorption because 
of intervening DLA systems. For example, the fluctuations in 
the number of X-ray photons received form different locations 
in 4C 41.17 exceeds roughly 100%. The main reason for this is 
small number statistics. The mean number of counts per pixel 
is ^ a few, and therefore Poisson noise will make variations be- 
tween individual pixels large. In 3C 294 (Fabian et al. 2003), 
the mean number of counts per pixel is higher by a factor of a 
few and therefore Poisson noise is less severe, although it is still 
large enough to introduce errors well in excess of 10%. Large 
holes in the X-ray map of 3C 294 are clearly seen, from which 
no X-rays at all are detected, and which are too large to arise 
from an intervening absorber. They probably arise from intrin- 
sic spatial variations in the population of relativistic particles 
responsible for the scattering. We next describe two methods 
to distinguish the intrinsic surface brightness fluctuations from 
those caused by absorption. 
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Fig. 3. — Power-law spectrum, Fe oc of an unobscured (dotted-red 
line), and an obscured X-ray source. The X-ray source is obscured by a col- 
umn of hydrogen corresponding to tq.s = 0.2 (dashed-hlue) and tq 5 = 0.3 
(solid-black). In Figure 1 the required column of hydrogen to produce 
To, 5 = 0. 1 can be read off as a function of z- The two mock data points show the 
error on the spectrum if such a source were observed and its spectrum binned 
in three energy bands: 0.3-0.6, 0.6-2.0 and 2.0-8.0 keV, assuming the the total 
number of photons received in the energy range 0.3-8.0 keV is 300 (see text). 



3.1. Spectral Imprints of DLA Systems 

Because the absorption cross section is a strong function of 
energy (crtot « E~^), the spectrum of an attenuated portion of the 
extended X-ray source can, in principle, reveal whether X-rays 
have been absorbed by an intervening column of gas. In order to 
make a crude estimate of the sensitivity of X-ray observations 
to DLA attenuation, we here quantify the absorption signature 
in a simple toy-model. In reality, an X-ray observatory has 
a potentially complex, energy dependent, detection efficiency, 
which must be taken into account in a more detailed analysis 
(more discussion on this is deferred to §3.1.1). 

In Figure 3, we show the spectrum of an unobscured X- 
ray source with a power law spectrum Fe oc with a = 1.3 
(dotted-red line , motivated by the inverse Compton spectrum 
observed in 4C 41.17). The units of the flux are arbitrary. Also 
shown is the spectrum of the same source, obscured by an inter- 
vening column of hydrogen with ro.5 = 0.2 (dashed-blue line) 
and T() 5 = 0.3 {solid black line). In these cases, absorption has 
removed ^ 18% and ^ 26% of the 0.5 keV photons, respec- 
tively. The two mock data points with errorbars are discussed in 
§3.1.1. Clearly, the spectra are most distinct at energies E < 0.6 
keV. 

Unfortunately, the unobscured case is purely hypothetical 
since omnipresent Galactic HI will always absorb a small amount 
of X-rays. The X-ray optical depth through the Galaxy rcai is 
typically 0.2, although there are large parts of the sky in which 
it is less than that by a factor of a few. Spatial variations in the 
total column of Galactic hydrogen across angular scales compa- 
rable to the sizes of known extended X-ray sources are negligi- 
ble. Therefore galactic HI will obscure the entire X-ray source 
by the same known amount and effectively steepen its spectrum 
at low energies (see §4.1). To find an absorption signature of a 
To 5 = 0.1 DLA system, for example, we therefore need to dis- 
tinguish between the rcai = 0.2 and rGai + To.5 = 0.3 cases (which 
are both shown in Fig. 3). 

In practice, for any source, one would proceed by first deter- 
mining the mean observed spectrum of the entire extended X- 
ray source. This mean spectrum can then be used as a template, 
to seek differential spectral deviations in spatial bins along the 
surface of the source. 

To quantify the difference in the spectra, we introduce the ra- 
tio TZ, which is the ratio of the number of photons received in 
the very soft (0.3-0.6 keV) and the soft (0.6-2.0 keV) band. In 
Figure 4, we plot TZ in the case of no absorption for a power 
law Fe oc spectrum as a function of a together with TZ 
for cases in which ro.5 = 0.2,0.3 and 0.4. Note that the case 
T0.5 = (dotted-red line) is hypothetical, but shown for com- 
pleteness. The case T0.5 = 0.2 (dashed-blue line) illustrates the 
case in which there is Galactic absorption only, whereas the 
cases ro.5 = 0.3 (solid-black line) and 0.4 (long-dashed-green 
line) correspond to cases in which there is additional absorption 
from a DLA system with T0.5 =0.1 and 0.2. 

Figure 4 shows, for example, that for a = 1.3 and T0.5 = 0, 
7^= 1.8 decreases to - 1.2/1.1 and 0.9 for T0.5 =0.2/0.3 and 
0.4, respectively. There is, of course, a degeneracy with a; for 
example, TZ can also equal 0.9 with no additional DLA absorp- 
tion when a = 0.9. Two approaches to lift this degeneracy are: 
(1) determining a via extrapolation from the £ > 1.0 part of 
the spectrum, and/or (2) computing 7?. as a function of position, 
which may reveal areas where TZ is significantly lower than in 
the rest of the source. The latter approach is similar to mapping 
out dust using an optical map, with the difference that obscured 
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parts of the X-ray source look 'bluer' , rather than redder. 
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Fig. 4. — Ratio of the number of photons in the observed 0.3-0.6 keV and 
0.6-2.0 keV bands, defined as H. as a function of the slope of the spectrum a. 
The different curves correspond to different amounts of absorption for a range 
of 0.5 keV optical depths of 0-0.4, as shown by the labels. The three mock 
data points denote the error on Ti and a for a = 0.8, 1 .3, and 1.8, provided the 
total number of detected photons in the energy range 0.3 — 8.0 keV is 300. 

3.1.1. Photon Statistics 

We next address the uncertainty in a and TZ given that we 
receive Wph photons in the energy range 0.3-8.0 keV from a spe- 
cific location on the sky. We bin the spectrum in the 0.6-2.0 and 
2.0-8.0 keV bands to two data points from which we extract the 
slope of the spectrum, a. The reason we only use these two 
bands is that they are not strongly affected by obscuration (see 
Fig. 3). 

As an example, we consider the case a= 1.3: the number 
of photons in the 0.3-0.6, 0.6-2.0 and 2.0-8.0 keV bands is 
149(Mph/300), 123(nph/300) and 28(nph/300), respectively for 
To 5 = 0.2, which represents the case in which there is Galactic 
absorption only. The ratio of the number of photons in the 0.6- 
2.0 and 2.0-8.0 keV bands is set by a. Assuming Poisson noise 
on the number of photons in each bin, we find the uncertainty 
in a and TZ to be, CTc ~ 0.2 and a-ji ^ 0.15. We repeated the 
calculation of the uncertainty for a = 0.8 and a= 1.8. The re- 
sults for a = 0.8, 1.3 and 1.8 are shown as three 'data points' 
in Figure 4. We find the following rough scaling relations of 
the errors: ^ 0.17(a/1.3) and ctk - 0.15(a/1.3)' l The 
a-dependence of the errors reflect that as a increases, the num- 
ber of photons in the harder bands decreases, which makes the 
uncertainty there larger. 

Figure 4 shows that the tq.s =0.3 and 0.4 curves are separated 
by 1 and 2a, respectively,from the tq.s = 0.2 curve, with a weak 
dependence on a. We write the level significance S at which 
absolution can be detected as: 

The above estimates ignore the processing of the input signal 
through the X-ray instrument. The typical response function of 



X-ray telescopes like Chandra are strong functions of energy. 
Chandra is a factor of ~ 1.5-2 and ~ 3 - 4 more sensitive at 
£ = 1.0 keV than it is at £ = 0.5 and 0.3 keV, respectively. In 
addition, Chandra's sensitivity decreases roughly linearly for 
E > 4.0 keV. Since the response function is known, one can cor- 
rect for this, leaving the values of TZ we calculated unchanged. 
The errors, however, are affected by this response function. For 
a fixed number of detected photons in the energy range 0.3- 
8.0 keV from a source with a certain a, we find that the total 
number of photons in the 0.3-0.6 and 2.0-8.0 keV bands is 
reduced by a factor of (0.7-0.8) and (0.8- 1.0), respec- 
tively , whereas it is increased by a factor of 1 .3 - 1 .4 in the 
0.6-2.0 keV band. We find that this barely affects the errors on 
TZ and a and conclude that equation (3) remains unchanged by 
a proper treatment of the energy response function. 

For comparison, Bechtold et al. (2001), detected 198 photons 
the 0.4-8.0 keV energy band for a source with a = 1.5 ±0.2 
in a 3 ks observation of Q1331H-170 with Chandra. First we 
notice that the error on a is the same as our estimate (note that 
198 photons in the energy band 0.4-8.0 keV, corresponds to 271 
photons in the band 0.3-8.0 keV for a = 1 .5). With this number 
of photons, their non-detection puts an upper limit of ^ 10^^ 
cm~^ to the column of hydrogen for the DLA system along the 
LOS at z = 1.77 at the 3a level. Assuming Z = O.IZ©, this 
can be translated to an upper limit to the value of ro.g for the 
DLA system of 0.26. If we interpret this upper limit of Tcai + 
To 5 = 0.46 at the 3 - cr level as a statement that the separation 
between the ro.g = 0.46 and 0.2 curves is 3a, then this is also in 
good agreement with equation (3), which predicts S = 2.3a. The 
numbers lie even closer together when we take into account that 
the value of rcai in the direction of Q1331H-170 was actually 
0.12 instead of 0.2, confirming that our simple estimate is a 
good first-order approximation for the expected signal. 

3.1.2. Required Integration Tunes on Chandra 

In §3. 1 . 1 we found that we need several hundred X-ray pho- 
tons to have a 2a detection of a foreground DLA system along 
the LOS, with the exact number depending on the slope of the 
X-ray spectrum of the background source, and the value of tq.s 
of the foreground DLA system. To make a map that spatially 
resolves a DLA system we need several hundred photons per 
pixel in the image. 

Ideally, we also need surface brightness maps with arcsecond 
resolution to be able to differentiate between compact clumps 
(<i < 5 kpc) and large disks {R > 15-30 kpc, or D > 30-60 
kpc, e.g. Chen 2005), since 1" corresponds to ^ 7-8.5 kpc in 
the redshift range 0.7-4.0. To reduce telescope time, the maps 
can be binned to a resolution of 2" x 2". If the DLA gas resides 
in a large (D > 30- 60 kpc) disk, it will be resolved in > 5 - 16 
resolution elements. 

To obtain > 3 x 10^ photons per resolution element, we need 
> 75 photons arcsec"^. We can express the required integration 
time for a detection in one resolution element (pixel) at signifi- 
cance level S as: 

^^t =3.0(^)'(f^)(|i^)(^)'Ms (4) 

= 3-0 (^) Ms, (5) 

where counts s"' is the count rate of the total number of pho- 
tons in the energy band 0.3-8.0 keV per arcsec"^. To arrive 
at the second fine (eq. 5), we assumed Z = O.IZ©, and used 
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equation (1) to eliminate T0.5. The strong redshift dependence 
of the required integration time puts constraints on which sub- 
classes of DLA systems can be detected in this manner. If 3 
Ms is taken as an absolute upper limit on the total integration 
time, the minimum required column of hydrogen in the DLA 
system as a function of redshift is A^m = 2.7 x 10^' [(l+z)/2]^'^ 
cwT^ in order for it to be detected. If we assume A^m cannot 
exceed 10^^ cm"^, then the method fails to detect DLA systems 
at z > 2.30. Note that although the maximum observed log(A^Hi 
/cm"^) = 21.7, theoretically there is no evidence for an upper 
limit at this value of A^hi- The numbers given above are all de- 
pendent on the assumed metallicity, on a and on the brightness 
of the background X-ray source. 

The quoted count rate of 2.5 x 10"^ counts s"' is motivated by 
the few recently observed high-redshift extended X-ray sources. 
For example, 4C 4 1 . 1 7 was observed for 1 00 ks on Chandra by 
Scharf et al. (2003) and has a few photons arcsec"^ over its 
^100 arcsec^ area. Another source, 3C 294, was observed for 
200 ksec with Chandra and has > 4-6 photons per arcsec^ over 
most of its area of ~ 200 arcsec^ (Fabian et al. 2003). Note that 
there are regions where the count rate is higher. 

Grego et al. (2004) observed a distant (z = 0.8) galaxy clus- 
ter for 100 ksec using Chandra. In their observation an area 
as large as ~ 10^ arcsec^ had more than 2.5 photons arcsec"^ 
The origin of the X-rays in this case is mainly thermal, and the 
analysis given above for a power-law inverse Compton spec- 
trum does not apply. However, a cluster X-ray spectrum in re- 
ality is more complex, with many heavy element emission lines 
superimposed on a multi-temperature thermal Bremsstrahlung 
spectrum (Raymond & Smith 1977). A power law with photon 
index ofawO.5 — 1 can therefore be used as a simple first- 
order approximation, and similar statistics will apply. An ad- 
vantage of using galaxy clusters (over the X-ray halos of radio 
galaxies) is their potentially larger angular size, which allows a 
search for the rarer DLA systems with more significant opacity. 

3.2. Combining X-ray and Sunyaev-Zeldovich Effect Maps 

Another way to determine whether a dimming of certain part(s) 
in the X-ray image is due to absorption or to an intrinsic bright- 
ness fluctuation in the source itself is to map the CMB temper- 
ature decrement, caused by the Sunyaev-Zeldovich (SZ) effect, 
along the surface of the source. This decrement should scale 
with the intrinsic X-ray brightness, regardless of whether the 
origin of the X-rays is thermal or non-thermal. In particular, 
thermal X-ray emission scales with the line-of-sight integral 
of Mg, whereas inverse Compton emission and the SZ surface 
brightness both scale with the LOS sight integral of n^ (where 
He is the electron number density). Local density variations can 
therefore cause thermal X-ray and SZ brightness variations. 
However, the thermal X-ray flux is unUkely to drop noticeably 
without an accompanying decrease in the SZ decrement. In 
contrast, an intervening DLA system would not have any effect 
on the SZ decrement. An SZ map with approximately arcsec- 
ond resolution would therefore be very useful in constraining 
the intrinsic X-ray brightness fluctuations. This resolution will 
be feasible with future millimter telescopes, such as the Ata- 
cama Large Millimeter Array (ALMA)^. Whether ALMA is 
sensitive enough to detect SZ temperature decrements due to 
sources like 4C 41.17 can be estimated as follows: As a first 
order estimate, the temperature decrement at ly ^ 100-200 
GHz can be approximated by ATcmb ~ (A^c 7^/10^*^ K cm"^) 

^See http://www.alma.nrao.edu/ 



/iK, where A',, is the column density of electrons in the X-ray 
source, Te = E/k is the electron temperature, and E is the mean 
electron energy. The product NeTg is directly proportional to 
the observed rate of X-ray photons, Nj, which is given by = 
edV/ (4-ndl(z)dilEj) photons cm"^ s"' sr"^ Here e = 5.4 x 
lO'^^TgNeil + z)'^ ergs s~' cm"^ is the volume emissivity from 
inverse Compton emission (e.g., Katz et al. 1996); diiz) is the 
luminosity distance; E^ is the mean observed photon energy, 
and dV and dil are differential volume and solid angle ele- 
ments, respectively. For 4C 41.17, ^ 2 photons arcsec"^ are 
observed, implying NeT,, ^ 10^*^ in the above units. Therefore, 
we find ATcMB ~ 10^ /xK, which can be detected by ALMA 
within several hours of integration (Kocsis et al. 2005). 

4. DISCUSSION 
4. 1 . The Absorber's Location along the LOS 

We found that it is possible to detect DLA systems in ab- 
sorption against background X-ray sources, but generally the 
redshift of the absorber along the LOS is not constrained (ex- 
cept in the rare cases in which the absorption signature is very 
strong; Bechtold et al. 2001). In this section we discuss possi- 
ble contaminants that could leave an absorption feature similar 
to a DLA system along a random position along the LOS. The 
sequence in which the contaminants are discussed is based on 
their distance to us. We start with Galactic HI and end with 
neutral gas physically associated with the source. 

Galactic HI is present in large columns; the average column 
of hydrogen within our own Galaxy is (AHisin|i»|) = 3 x lO^'' 
cm"^, where b is Galactic latitude and (...) denotes the average 
over an annulus at b Lockman e.g., 2003. Assuming a solar 
abundance pattern of Z = Q.IZq (which is the appropriate value 
for the interstellar medium (ISM), e.g. Wilms et al. 2000) for 
this column, this translates to (TGaisin|^|) = 0.19. Because of 
spatial variations of Galactic A^hi. there are large portions of the 
sky where roai is lower by a factor of a few. These spatial vari- 
ations occur predominantly on large (> 1°) angular scales, and 
therefore rcai can be considered as a constant over the entire X- 
ray source. The rms spatial fractional variation in Galactic A^hi 
is < 3% on angular scales < 1°, with larger fluctuations occur- 
ring on smaller angular scales only in directions of anomalously 
high Galactic A^hi (Lockman 2003). This leads to the conclu- 
sion that the Galactic column can be treated as constant across 
the extended X-ray source, and is unUkely to present a problem 
beyond changing the effective slope of the X-ray source at low 
(£■ < 0.6 keV) energies. 

The regions of anomalously high A^hi are not found only in 
the Galactic plane, but clusters of HI clouds associated with the 
galaxy are found at Galactic latitudes > 10-20° (Lockman & 
Pidopryhora 2004). The association of this gas with the galaxy 
is clear because its kinematics closely matches that of the planar 
gas. Less clear is this association for the High Velocity Clouds, 
which have velocities that deviate by > 50 km/s (Wakker & 
van Woerden 1997) from the aUowed range based on a model 
of Galactic rotation. The nature of these HVCs is not well un- 
derstood, mainly because the distance to them is observation- 
ally poorly constrained. Possible contamination by HVCs or 
extraplanar galactic clouds can be avoided simply by avoiding 
the directions in which they are detected, as mentioned earlier 
However, potentially smaller versions of HVCs or any kind of 
extraplanar HI cloud may contaminate our signal. 

A smaU compact HI cloud with an angular size of ~ 5" and 
a column of a few times 10^" cm~^ would go undetected in cur- 
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rent HI observations, and could mimic the absorption signature 
of a DLA system at a cosmological distance. We discuss be- 
low why these clouds, if they somehow manage to form, are 
probably short-lived and probably not a contaminant we should 
worry about. For a cloud with A^hi = 3 x 10^" cm"^ at a distance 
Z) = 100 kpc with an angular size of ^? = 5 arcsec, we obtain 
an HI mass of ~ 10(D/100 kpc)^ (0/5")^ compressed in a 
cloud of radius 1.2 (Z)/100 kpc) {9/5") pc. With a mean num- 
ber density of ~ 2.5 x 10^(^/10 Mpc)"' (6'/5"r' cm'^ these 
clouds are only self-gravitating if the temperature inside them 
is ~ 4(D/100 kpc) K. These cold, high density environments 
are typical for molecular clouds ( although these are more mas- 
sive, with M ^ IO^Mq). Therefore one would expect these 
clouds to be sites of starformation as well, in which one cloud 
would form only a few stars (or one star) on a short timescale. 
Note that if these clouds were truly isolated, another threat to 
their stability would be thermal contact with a hot surrounding 
medium (the Galactic halo gas or the intergalactic medium). 
Heat conduction would rapidly raise the cloud's temperature, 
which would cause it to expand and lower it's column density. 

The required columns are also typically reached in cold disks 
of late-type galaxies in the nearby universe, which would be 
bright and easy to identify in optical images. Furthermore, even 
if the HI is associated with faint galaxies, it will be easier to 
identify than in the case of DLAs discovered in optical quasar 
spectra, because there would be no bright background quasar 
dominating the optical images. 

Another possibility is that the hydrogen gas is physically as- 
sociated with the source. Whether a significant amount of neu- 
tral hydrogen can co-exist with the hot. X-ray emitting gas 
should depend primarily on the local cooling efficiency. In re- 
cent work, Birnboim & Dekel (2003) and Keres et al. (2004) 
have argued that significant neutral gas is present only around 
low-mass, high-redshift galaxies. Sources that are sufficiently 
X-ray bright to allow the detection of a DLA system are Ukely 
to be associated with more massive objects (such as powerful 
radio galaxies), in which the local neutral hydrogen would be 
absent. In any case, absorption in a region that is several arc- 
seconds across would imply that the neutral hydrogen gas asso- 
ciated with the source itself would extend several tens of kilo- 
parsecs across, which would be an interesting alternative inter- 
pretation of any future detection. 

We conclude that the two contaminants that cannot be ruled 
out observationally, the compact, cold. Galactic clouds or mini 
HVCs and the HI associated with the X-ray source, both offer 
interesting alternative interpretations of any future detection. 

4.2. The Uncertainties in fiNm,z) 

As mentioned above, an important uncertainty in our predic- 
tions is the poorly known column-density distribution /(A^m, z) 
at low redshifts and high columns. The total number of DLA 
systems used by P03 to constrain /(A^hi,z) in the DLA range of 
Nm was 1 14. These data were binned into four redshifts bins 
covering the redshift ranges z = 0.0-2.0, 2.0-2.7, 2.7-3.5 
and > 3.5 each containing between 23 and 49 DLA systems. 
P03 fitted a Schechter function, /(Nm , z) = (/* /N^) (Nm/N^T^ 
exp(-A^Hi/A^*), to the data in each bin. Prochaska & Herbert- 
Fort (2004) identified six new DLA systems at z > 3 .5 in an au- 
tomated search for DLA systems in the quasar spectra of Data 
Release 1 of the Sloan Digital Sky Survey (SDSS), and found 
that including these six new DLA systems changed the fitting 
parameters (/* , A^* , /?) of P03 significantly at z > 3 .5 which they 



attribute to the small sample size in this bin^. For the same rea- 
son the obtained values for the fitting parameters in the lowest 
redshift bin are uncertain. 

Not only the fitting parameters, but also the fitting function it- 
self is uncertain. Using a Schechter function to describe /(Nm , z) 
introduces a sharp cutoff at high A^hi- P03's /(A'hiiZ) provides 
a good fit to the colunm density distribution over most of the 
z and A'hi range. However, at the high-A'Hi regime that is rele- 
vant for our purposes, other functional forms (with less sharp 
cut-offs in A^Hi) can provide equally good fits. 

We therefore find that the largest uncertainties in /(A'hi , z) are 
in regions of (A'Hi,z)-space that are most relevant for our pur- 
poses. As can be seen in Figure 1 and equation (1), the column 
required to produce a certain T0.5 increases rapidly with redshift 
(oc (1 +z)^'^, for Z = O.IZ0). Because the high column systems 
are rare, the largest contribution to /cov in the P03 model comes 
from low redshift. As mentioned in §2.4, 75% of the contribu- 
tion to /cov comes from z < 0.6. 

This dominance of low redshift DLA systems is artificial 
since the current f(Nui,z) does not capture any evolution in the 
number density of DLA systems in the redshift range z = 0-2. 
Observationally, it is known that the number density of DLA 
systems increases with redshift between z = and z = 1 (Rao & 
Turnshek 2000). An improved determination of fiNHi,z) that 
does capture the redshift evolution in the range z = - 2 will 
weaken the dominance of low redshift systems. 

In summary, our estimate of /cov is uncertain, but is likely to 
be conservative, since P03 may have underestimated the num- 
ber of high A^Hi systems, especially in the highest redshift bin. 
In addition , the degree to which low redshift DLA systems 
dominate the contribution to /cov is exaggerated with the cur- 
rent /(A'hi, z). 

4.3. Using X-ray Point Sources to constrain /(NuItZ). 

An intriguing possibihty is to exploit the enhanced sensitivity 
to low-redshift systems to constrain /(A^hi,z) at low redshifts. 
The brightest X-ray sources Bechtold et al. (2001) have count 
rates of ~ 1 photon s~' ; ~ 3 - 1 5 ks integrations on Chandra can 
be used to constrain whether and how much additional (to the 
Galaxy) absorption is present in the spectra of these sources. 
X-ray sources that are fainter by factor of ~ 10, which are more 
common, can detect a foreground DLA system with T0.5 = 0.1 
at the 3(7 level in 25 ksec according to equation (5) (Note that 
for point sources we need to multiply eq. [5] by 4, since in its 
derivation we assumed the X-ray was extended and binned into 
2" X 2" pixels.) The number of X-ray sources of this brightness 
or more is expected to be a few tens per square degree (e.g. 
Moretti et al. 2003). Although the redshift of the absorber is not 
known, with the exception of the cases in which the absorption 
signal is very strong (Bechtold et al. 2001), it is constrained by 
the redshift of the source. 

It therefore seems promising to look at archival data in which 
sources brighter than 0.1 counts s~' have been observed for 
25 ks, or more generally, any archival data that has > 2.5 x 
10^(5/3(t)^ counts from a given location and see if any addi- 
tional absorption is detected at the level S. The main difficulty 
with this method is that the X-ray point source itself may be 
obscured. As a result, the amount of absorption attributed to 
any DLA system will be uncertain. Statistically, the distribu- 
tion /(A'hi 7 z) of DLA systems could still be inferred, but would 

' Prochaska & Herbert-Fort (2004) also note that P03 have obtained systemat- 
ically too low values of Nm for logCiVm/cm"^) > 21. 



8 



require modeling of the intrinsic obscuration (e.g. Gilli et al. 
(2001)). A more robust approach would be to select unobscured 
sources, and derive upper limits on any DLA absorption, which 
would constrain (put upper hmits) on f(Nm,z). 

4.4. Constraints on Galaxy Formation Scenarios 

As mentioned in § 1, a spatially resolved map of the DLA gas 
could distinguish between the two models in which the gas re- 
sides either in a large disk, or in smaller cold clumps accreting 
onto to the galaxy. The large-disk models apply at z > 2, and as 
we saw in §2.3 the contribution to /cov form z > 1.4 is negligi- 
ble. Although we found in §4.2 that this dominance of the lower 
redshift DLA systems is exaggerated, we also found in §3.1.2 
that long, > 3 Ms integration times are required to detect DLA 
systems at z > 2.3 (eq. 5). Unless future X-ray observations 
reveal extended X-ray sources with count rates > 2.5 x 10"^ 
counts s"\ we will not be able to detect z > 2.3 DLA systems 
in absorption with log(A^Hi/cm""^) ^ 22, which would exclude 
all currently known DLA systems. 

We conclude that current X-ray telescopes most likely will 
not be able to produce maps that would enable us to distin- 
guish between the distributions of gas predicted by the dif- 
ferent galaxy formation scenarios. This is even more so, be- 
cause the integration times mentioned above are for X-ray maps 
smoothed to 2" x 2". At z = 2, 2" corresponds to ~ 16 kpc, and 
small clumps would not be resolved. It is, however, possible to 
produce maps of DLA gas at z < 2. Knowledge of the z < 1 sys- 
tems is limited, because currently only 23 are known. For these 
23 DLA systems, 14 have confirmed galaxies as counterparts, 
which are consistent with being drawn from a field population 
of galaxies (e.g. Chen 2005). For a few of these galaxies with 
a stellar disk, the rotation curve of the ionized gas within the 
stellar part of the galaxy was compared with the velocity differ- 
ence between the galaxy and the DLA gas. Based on the results 
and the distance of the DLA gas to the center of the galaxy, it 
was found that the gas in the z < 1 systems was consistent with 
being in a large disk, extending to radii larger by a factor of \/2 
than those observed in local disk galaxies (see Chen 2005 and 
references therein). Producing maps of the DLA gas in these 
z < 1 systems may provide an important test for this result. In- 
creased knowledge of the z < 2 systems will clearly add to the 
understanding of their higher redshift counterparts. 

5. MAPPING OUT KNOWN DLA SYSTEMS 

So far we have focused on examining an extended X-ray 
source and searching for evidence of an intervening DLA sys- 
tem. As mentioned above, several hundred DLA systems are 
known to exist in the redshift range z = 0.1-4.6. A complemen- 
tary approach would be to look for extended X-ray sources be- 
hind these DLA systems, and examine these extended sources 
for evidence of extended absorption from the DLA system. 

Scharf et al. (2003) choose 4C 41.17 for X-ray study be- 
cause it was an example of a cluster of submillimeter sources, 
centered on a powerful radio galaxy, and was therefore consid- 
ered a probable proto-galaxy-cluster with X-ray emission. Ex- 
tended X-ray sources at high redshifts may indeed turn out to 
be common, especially around radio-loud galaxies. Of the sev- 
eral hundred known DLA systems in front of quasars, ^ 10% 
w be along sight-hnes to radio loud quasars. A certain fraction 
of these radio loud quasars will be similar to 4C 41.17. This 
fraction is unknown at present, but is unlikely to be negligible. 



The first question is, what fraction of the known DLA sys- 
tems have To.5 > 0.1, and would therefore be good candidates 
for such a study. A table containing ~ 300 DLA systems is 
given in Curran et al. (2002), for which we calculated the cor- 
responding values of tq s- Since the metallicity is unknown 
for most of them, we assumed it to be O.IZq. We found a 
large spread in the values of tq.s, with the maximum as high 
as To.5 = 1-2. We also found that 12 and 7 DLA systems have 
To.5 > 0. 10 and > 0.20, respectively (while this number was re- 
duced to 8/6 when we lowered the metallicity to O.OIZq). The 
DLA systems with tq.s > 0.20 all he at z < 0.7, and there are 
three at z > 1.4 with tq.s = 0.10-0.12. The highest value of 
ro.5 = 1.2 is for the source detected by Bechtold et al. (2001), 
already mentioned above. The spectral signature is so evident 
in this case, that, based on the amount of X-ray absorption and 
the known HI column in the DLA system from its radio absorp- 
tion signature, the metallicity of the DLA gas could be derived, 
and was found to be ^ 0.23Zq. 

The second question is: what fraction of these relatively opaque 
absorbers have extended X-ray emission surrounding the opti- 
cal quasars in which they were discovered? An indicator for 
the presence of extended X-ray emission may be radio loud- 
ness (powerful radio galaxies may drive such X-ray emission, 
such as that found in the case of 4C 41.17). The table of DLA 
systems composed by Curran et al. (2002) lists the radio flux 
densities of the background quasars. We find that among the 12 
DLA systems with T0.5 > 0. 1 , at least three are in front of either 
a 3C or 4C source, which are radio loud. Among the other 8, 
at least 6 have detections in the radio. An intriguing example 
is Q1354-H258, which is similar to 4C 41.17: it appears to have 
an extended structure in the radio, and is additionally known to 
have a diffuse Lya blob around it (Heckman et al. 1991). This 
source would be a good first candidate to look for accompany- 
ing extended X-rays. The DLA system in front of Q1354H-258 
is located at z = 1-42, and it has a column of 3.2 x 10^' cm~^ 
(e.g., Rao & Turnshek 2000). Note that the metallicity of this 
DLA system is Z - O.O2Z0 ( Vladilo 2002) instead of Z = 0. IZq 
which slightly lowers ro.5 from 0.12 to 0.09. This implies that 
if extended X-rays are detected from Q1354-I-258, it should be 
brighter than 4C 41.17 for a detection at the 5 = 2-(t level in 
^ 3 Ms on Chandra. 

The above exploratory remarks suggest that a full systematic 
search for DLA systems, whose parent quasars are similar to 
4C 41.17, and deep X-ray observations of these sources, is a 
promising alternative approach to discover X-ray silhouettes 
of DLA systems. 

6. CONCLUSIONS 

The neutral hydrogen column densities in observed DLA sys- 
tems are in the range log(A?Hi/cm~^) = 20-21.7 (e.g. Cur- 
ran et al. 2002) and are high enough to absorb a measurable 
fraction of soft X-ray photons from a background source. We 
defined the parameter T0.5, which is the optical depth through 
DLA systems for a photon with an observed energy of 0.5 keV. 
We showed in §2.3 that ~ 1% of the sky is covered by DLA sys- 
tems which have ro.5 > 0.2 when the background X-ray sources 
lie at redshifts z 1.4. This result depends on the assumed col- 
umn density distribution /(NhijZ), which is at present poorly 
constrained, and also on the mean metalUcity of the DLA sys- 
tems, which we assumed to be Z = 0. IZ0 based on the observed 
value (e.g. Prochaska et al. 2003). We discussed uncertainties 
in /cov in §2.4 and found that the values we found are most 
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likely conservative. 

An optical depth of T0.5 = 0.2 will introduce a dimming in 
the soft X-ray bands by ^ 20%, an amount that is much lower 
than intrinsic brightness fluctuations. However, as we discussed 
in §3, absorption by an intervening DLA system can be distin- 
guished from intrinsic brightness fluctuations by looking for the 
spectral signatures of absorption (§3.1 , Fig. 3) and/or combin- 
ing the X-ray observations with an SZ decrement map of the 
area (§3.2). We showed in§3.1.1thatin order to reli ably detect 
a DLA absorption signature on top of the spectral imprint left 
by Galactic HI, > 300 photons per angular resolution element 
are required in the 0.3-8 keV band to detect a DLA system with 
To, 5 = 0.2 at the 2cr level. 

For diffuse, high redshift X-ray sources similar to 4C 41.17 
(Scharf et al. 2003) and 3C 294 (Fabian et al. 2003), this re- 
quires observation times of a few megaseconds on Chandra 
when the maps are smoothed to ~ 2" resolution. The influence 
of the brightness of the background X-ray source, the slope of 
the X-ray spectrum and the value of ro.5 on the required inte- 
gration time on Chandra are detailed in §3.1.2. 

We discuss the implications of the poor constraints on the ab- 
sorber's redshift using this method in §4.1 and find two possible 
contaminants that cannot be ruled out observationally. These 
are (1) neutral hydrogen physically associated with the X-ray 
source and (2) a small, compact nearby HI cloud. Both are as- 
trophysically interesting as we argued in §4.1. We suggested 
a new method to use archival data of bright (count rates > 0.1 
photons s"') X-ray point sources to constrain the poorly known 
low-redshift column density distribution function, fiNm,z), in 
§4.3 

The inverse method is discussed in §5 in which we presented 
the possibilities of mapping out known DLA systems using an 
X-ray background source. We found that out of the sample of 

300 DLA systems tabulated in Curran et al. (2002), 12 have 
T0.5 > 0.1, with values anywhere in the range tq.s = 0.1-1.2. 
We found that at least 4 of these sources resemble 4C 41.17 
and 3C 294 in their radio loudness and at least one of them 
also has a surrounding extended Lya halo. These 4 sources are 
promising candidates to have extended X-ray emission associ- 
ated with them, and, if such X-ray emission is revealed in deep 
X-ray follow-ups, they would be suitable to search for DLA 
silhouettes. 

Although challenging. X-ray silhouettes are likely to be the 
only viable method to directly spatially resolve hydrogen gas in 
distant DLA systems in the foreseeable future. Among planned 
future radio telescopes, only SKA will be able to detect 21cm 
emission from neutral hydrogen at z > 0.5 (e.g., van der Hulst 
et al. 2004), and extended Lya emission would only be de- 
tectable for DLAs in the immediate vicinity of a luminous ion- 
izing source. 

The planned X-ray observatories, the X-ray Evolving Uni- 
verse Spectrometer (XEUS; Parmar et al. 2003) and Generation 
X (Zhang et al. 2002) are expected to have ~ 0.1" - 1" resolu- 
tion, and a sensitivity of 1 or 2 orders of magnitude better than 
Chandra, extending to energies as low as 0.1 keV. At these low 
energies, absorption by a foreground DLA system will leave a 
much stronger imprint. The X-ray silhouettes of DLA systems 
should be feasible to detect in routine observations with these 
planned instruments. 
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